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Abstract. 

This is a draft describing the calculation of the evolution of the gluon density in x from an initial 
value X = xq = 0.01 to smaller values, up to x = 10“® in the hard pomeron formalism with a running 
coupling introduced on the basis of the bootstrap equation. The obtained gluon density is used to 
calculate the singlet part of the proton structure function. Comparison with experiment and the 
results following from the fixed coupling evolution is made. 
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1 Introduction. 

Recent results obtained at HERA |Q, ^ may be interpreted as a manifestation of the hard pomeron, 
which naturally explains a sharp rise of F 2 {x,Q^) at low x. The original BFKL hard pomeron, 
however, has a drawback of treating the coupling constant as fixed, since it sums only powers of 
log 1 /x and not those of log . A rigorous way to introduce a running coupling into it still remains 
beyond the possibilities of the theory, since it inevitably involves a problem of low behaviour and 
thus of confinement. In a series of papers [|[ |[ ||, ^ we have adopted a more intuitive way to attack 
this problem, based on the so-called bootstrap relation IJ , which is, in fact, the unitarity condition 
for the t-channel with a colour quantum number of a gluon. It ensures that the one-reggeized-gluon 
exchange supposed to give a dominant contribution in this channel is unitary by itself, which is a 
necessary requirement to use it as an input for the construction of the BFKL pomeron. 

Assuming that this fundamental requirement should be preseved in the theory with a running 
coupling, we proposed a minimal modification of the BFKL pomeron equation which, on the one 
hand, satisfies the bootstrap condition and, on the other hand, leads to the standard results with the 
running coupling in the double log (DL) limit, i.e. when leading terms in the product log 1/a; log 
are summed. This modification reduces to the substitution of every momentum squared in the 
pomeron equation by a function r]{k) which at large k behaves as k"^/2as{k‘^). The behaviour of rjik) 
at low k remains beyond any theoretical controle. We parametrize it as 

r]{k) = (6o/87r)(fc^ -f m^) In ^ (1) 

where A is the standard QCD parameter, 6o = 11 — {2/3)Nf and the effective gluon mass m simulates 
both the confinement effects and the freezing of the coupling. 

Solving numerically the pomeron equation in this approach we found two supercritical pomerons 
Q. Adjusting the mass m to fit the experimental slope of the leading pomeron of 0.25 {GeVlc)~‘^ we 
obtained for their intercepts 

Ao = 0.384, Ai= 0.191 

and the slope of the subdominant pomeron results a'l = 0.124 {GeV/c)~2. Calculating observable 
quantities with only these two asymptotic states taken into account we found that the picture which 
emerges, in all probability, corresponds to energies much higher than the present ones |^. In particular 
the average {k±) was found to be very large (~ 10 GeV/c) and independent of energy, which may 
indicate a saturation of its growth observed at present energies [|| . 

To describe the present experimental data it is then necessary to take into account all the states 
from the pomeron equation spectrum. This can be achieved by converting the pomeron equation into 
an evolution equation in I /x and solving it with an initial condition at some (presumably small) value 
X = xq. In such an approach, taking a nonperturbative input at a; = a;o adjusted to the experimental 
data, also the problem of coupling the pomeron to the hadronic target is solved in an effective way. 

This note is devoted to realizing such a program. In Sec. 2 we state our basic equations. The most 
difficult part of the program is to pass from the gluon density to the observable structure function. It 
is discussed in Sec. 3. Sec. 4 is devoted to fixing the initial gluon distribution for the future evolution. 
In Sec. 4 we present our numerical results. Sec. 5 contains a discussion and some conclusions. 



2 Basic equations 

For the forward scattering amplitude the pomeron equation reads 
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{H - E)ip = ipo 


(2) 


Here ip is a semi-amputated (one leg amputated only) pomeron wave function; E = 1 — j is the 
pomeron ’’energy”, related to its complex angular momentun j] H = T + V is the ’’Hamiltonian” 
consisting of the kinetic energy given by the sum of the two gluon Regge trajectories, T = —2lo and 
of the potential energy V. With a running coupling introduced according to ii both are expressed 
via the mentioned function rj (Eq. (1)); 


and 


Vip{k) = - 


T{k) = 


T1T2 


fVe f d^k'-qjk) 

( 27 r)^ J r]{k')r]{k — k') 


drk'^l;{k') 
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(3) 


(4) 


{2'kY J \r]{k — k') r]{k)r]{k') 

where Nc is the number of colours and Ii( 2 ) are the colour operators for the two interacting gluons; 
in the vacuum channel we have T 1 T 2 = —Nc- Finally, the inhomogeneous term ipQ represents the 
interaction vertex between the pomeron and the hadronic target. 

Taking the Mellin transformation of (2) one converts it into an evolution equation in 1/x: 

d 


which should be supplemented with an initial condition at some x = xq 


(5) 


ip{xo,k) = ipo{k) 


( 6 ) 


containing the nonperturbative input about the coupling to the hadronic target. 

The physical interpretation of the pomeron wave function is provided by the fact that in the DL 
approximation Eq. (5) reduces to an equation for the fully amputated function (p{x, k) = rj{k)ip{x, k): 
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<P{x, k) 


3as(fc^) 
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(p{x,k) 


(7) 


which coincides with the standard equation for the unintegrated gluon density xg{x, k^) in the DL 
limit. In fact, this circumstance lies at the root of our method to introduce a running coupling into 
the scheme. Thus we may identify 

(p{x, k) = cxg{x, k) (8) 


The normalizing factor c cannot be determined from the asymptotic equation (7). We shall be able 
to fix it by studying the coupling of the pomeron to the incoming virtual photon in the next section. 


3 Coupling to the virtual photon 

Once the function (p proportional to the gluon density is determined, one has to couple it to the 
projectile particle to calculate observable quantities. In particular, to find the structure function of 
the target one has to couple the gluons to the incoming virtual photon, that is, to find the colour 
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density p{q, k) which connects the photon of momentum q to the gluon of momentum k. This problem 
is trivial within the BFKL approach with a hxed small coupling. Then it is sufficient to take the colour 
density in the lowest order, which corresponds to taking for it the contribution of a pure quark loop 
into which the incoming photon goes po{q, k). 

The problem complicates enormously when one tries to introduce a running coupling into p. Then 
one has to take into account all additional gluon and qq pair emissions which supply powers of the 
logarithms of transverse momenta. Apart from making the coupling run, they will evidently change 
the form of p{q, k). Unfortunately the bootstrap relation can tell us nothing about the ultimate form 
of the colour density with a running coupling, which essentially belongs to the t-channel with a vacuum 
colour quantum number. So we have to find a different way to introduce a running coupling into p. 

A possible systematic way to do this consists in applying to the photon-gluon coupling the DGLAP 
evolution equation. One may separate the colour density from the rest of the amplitude by restricting 
its rapidity range to some maximal rapidity j/o ~ log Q'^ (which, of course should be much smaller than 
the overall rapidity Y ^ logQ^/x). Then the kinematical region of p{q,k) will admit the standard 
DGLAP evolution in Q^. Solving this equation one will find the quark density at scale of the gluon 
with momentum k (i.e. essentially the structure function of the gluon with the virtuality fc^). This 
is exactly the quantity needed to transform the calculated gluon density created by the target into 
the observable structure function of the target. As a starting point for the evolution one may take 
the perturbative colour density po at some low when the logs of the transverse momenta might be 
thought to be unimportant. 

This ambitious program, combining both evolution in both 1/x and Q^, does not, however, look 
very simple to realize. As a first step, to clearly see the effects of the introduction of a running coupling 
according to [ 3,4], we adopt a more phenomenological approach here, trying to guess a possible correct 
form for p{q, k) on the basis of simple physical reasoning and also using the DL approximation to hx 
its final form. 

With a pure perturbative photon colour density one would obtain for the 7 *p cross-section 



(Pkpo{q, k)4>{x, k) 
{2TT)‘^rf'{k) 


(9) 


In fact, the projectile particle should be coupled to the full pomeronic wave function (ji/rj^■ From the 
physical point of view this expression is fully satisfactory for physical particles. However it is not for 
a highly virtual projectile. 

To see this, we first note that for the forward amplitude our method of introducing a running 
coupling reduces to a very simple rule: the scale at which the coupling should be taken is given by 
the momentum of the emitted real gluon {{k — k'Y va the upper rung in Fig. 1). Now take very 
large and apply the DL approximation. Then the momenta in the ladder become ordered from top to 
bottom 

>> k"^ » k''^ » . 


In this configuration, as can be traced from (2) and (9), all a^’s acquire the right scale (i.e. corre¬ 
sponding to the DGLAP equation) except for the upper rung: as{k‘^) appears twice. This defect can 
be understood if one notices that the upper gluon is, in fact, coupled to a virtual particle. If this 
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particle were a gluon, then the interaction (4) would cancel one of the two a(fc^)’s and substitute it 
by an a taken at the scale corresponding to its own virtuality. We assume that something similar 
should take place also for virtual quarks to which the gluon chain may couple. The scale of the particle 
momenta squared which enter the upper blob in Fig. 1 should have the order (this is the only scale 
that remains after these momenta are integrated out). As a result the lowest order density should be 
rescaled according to 


Po{q,k) 


<Ql) 


Po{q,k) 


( 10 ) 


a{k‘^) ‘ 

where has the same order as Q^. 

The approximation we assume in this paper is that the substitution (10) is sufficient to correctly 
represent the photon colour density with a running coupling. We shall check its validity by studying 
the quark density which results from (10) in the DL approximation and comparing it with the known 
result based on the DGLAP equation. 

Explicitly the zeroth order density po has the following forms for the transverse (T) and longitudinal 
(L) photons (see e.g. Q and do the integration in the quark loop momenta) 


(a^ + (l- 

afm + 2z^)g{z) - I) + ^ . Al - g{z))) 

a{l-a) + C ) 

(11) 

3e^ 

f 

Zf / da 9 {z)) 

^ Jo Cl(l ~ Ck) + C 

(12) 


Here the summation goes over the quark flavours. The dimensionless variables C and z are defined as 


t f 

^ = OS’’ ^ “ 


1 


Q2 ’ - 4Q2 a{l -a) + C 

and TO/ and Zf are the mass and charge of the quark of flavour /. The function g{z) is given by 


(13) 


aiz) = 


1 


^z‘^ + l + z 


2zVz^l^ Vz^ + l-z 
The structure function is obtained from the cross-section by the standard relation 

F2{x,Q^) = 

Tre^ 


(14) 


(15) 


In the DL limit only the transverse cross-section contributes. We can also neglect the quark masses 
in this approximation. Then, with a substitution (10), from (9), (11) and (15) we obtain an expression 
for the quark (sea) density of the target 

W / t-lnP I ■;«(« + (I - ») )((l + 2^ )9(») - I) (16) 

where g{z) is given by Eq. (14) and we assumed that large values of k'^ < contribute in accordance 
with the DL approximation. In this approximation the asymptotics of the gluon density xg{x, k^) and 
consequently of ipix, fc^) is known: 


(f>{x, k^) = cxg{x, k^) ~ cexp 


a In — In In k"^ 
X 


(17) 
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where a = 48/6o- Putting (17) into (16), after simple calculations we find the asymptotical expression 
for the quark density 

a In —Inin (18) 

X 

On the other hand, from the DGLAP equation we find, with the same normalization 

xq{x,k^) ~ ^ 

As we observe the approximation (10) for the colour density of the photon projectile leads to the 
correct relation between the quark and gluon densities in the DL limit. This justifies the use of (10), 
at least for high Xjx and . Comparing (18) and (19) we also obtain the normalization factor c 
which relates the pomeron wave function to the gluon density 

c = 7r^6o/3 (20) 


In In 
a In 1 /: 


■ exp 


a In — In In 

X 


(19) 


xq{x^ k^) ~ 


4c / In In k'^ 


7r^6n V a In 1/x 


exp 


4 The initial distribution 

To start the evolution in l/x we have to fix the initial gluon density at some small value x = Xq. 
Evidently, the smaller is xq, the smaller is the region where we can compare our predictions with the 
experimental data. On the other hand, if xg is not small enough, application of the asymptotic hard 
pomeron theory becomes questionable. Guided by these considerations we choose xg = 0.01 as our 
basic initial x although we also tried x = 0.001 to see the influence of possible subasymptotic effects. 

The initial wave function 4>{xg, k'^) has to be chosen in accordance with the existing data x = xg 
and all k'^ available. The experimental F 2 is a sum of the singlet and nonsinglet parts, the latter 
giving a relatively small contribution at a: = 0.01. Our theory can give predictions only for the singlet 
part (and one of the criteria for its applicability is precisely the relative smallness of the nonsinglet 
contribution). The existing experimental data at a; = 0.01 give values for F 2 averaged over rather 
large intervals of x and Q^. For all these reasons, rather than to try to adjust our initial (j){xg,k‘^) 
to the pure experimental data, we have preferred to match it with the theoretical predictions for the 
gluon density and the singlet part of F 2 given by some standard parametrization fitted to the observed 
F 2 in a wide interval of and small x. As such we have taken the GRV LO parametrization [jl^ . 
The choice of LO has been dictated by its comparative simplicity and the fact that at a; = 0.01 the 
difference betwen LO and NLO is insignificant. 

Thus, for the initial distribution we have taken the GRV LO gluon density with an appropriate 
scaling factor. Putting this density into Eqs. (8),(9) and (15) one should be able to reproduce the 
sea quark density and thus the singlet part of the structure function. In the GRV scheme the relation 
between the gluon density and the quark density is much more complicated and realized through the 
DGLAP evolution. Since the DGLAP evolution and the pomeron theory are not identical, one should 
not expect that our initial gluon density should exactly coincide with the GRV one to give the same 
singlet structure funcction. One has also to have in mind the approximate character of our colour 
density p at small Q^. In fact, with the initial (j) given by (8) and the gluon density exactly taken 
from the GRV parametrization a,t x = 0.01 we obtain a 30% smaller values for the singlet part of the 









7 


structure function as given by the same GRV parametrization, the difference growing at low Q^. To 
make the description better we used a certain arbitrariness in the scale Q\ which enters (10) and also 
the scale at which the coupling freezes in the density p. The optimal choice to fit the low data is 
to take 

^ ^ ln((0.17 * Q 2 + 0.055 (Gey/c)2)/A2) 

With this a{Ql) the obtained singlet structure function at x = 0.01 has practically the same 
dependence as the GRV one, although it results 30% smaller in magnitude. This mismatch can be 
interpreted in two different ways. Either we may believe that the gluon density given by the GRV 
is the correct one and the deficiency in the singlet part of the structure function is caused by our 
approximate form of the colour density p (which is most probable). Or we may think that the colour 
density to be used in the DGLAP should coincide with ours only for large enough and 1/x and 
at finite values they may somewhat differ (our relation (8) was established strictly speaking only in 
the DL limit). Gorrespondingly we may either take the relation (8) as it stands and use the GRV 
LO gluon density at x = 0.01 in it, or introduce a correcting scaling factor 1.3 which brings the 
structure function calculated with the help of (9)-(15) into agreement with the GRV predictions. In 
the following we adopt the second alternative, that is we assume that our initial gluon distribution 
at X = 0.01 is 30% higher that the one given by the GRV parametrization. The singlet part of the 
structure function at x = 0.01 calculated from (9)-(15) with this choice is shown in Fig. 2 together 
with the GRV predictions. However one can easily pass to the first alternative by simply reducing our 
results by factor 1.3. 

5 Evolution: numerical results 

With the initial wave function cj){x = 0.01, k^) chosen as indicated in the preceding section we solved 
the evolution equation for 10“® < x < 10“^. The adopted calculational scheme was to diagonalize 
the Hamiltonian in (2), reduced to one dimension in the transverse momentum space after angular 
averaging, and represent the initial wave function as a superposition of its eigenvectors. To discretize 
k'^ a grid was introduced, after which the problem is reduced to a standard matrix one. To check 
the validity of the obtained results we have also repeated the evolution using a Runge-Kutta method, 
resulting in a very good agreement. The final results obtained for the gluon distribution xg(x, Q^) as 
a function of for various x are shown in Figs. 3 and 4 and as a function of x for various in 
Figs. 5 and 6. Figs. 3 and 5 correspond to x and presently available, whereas Figs. 4 and 6 show 
the behaviour of the calculated gluon density in the region up to very small x and very high well 
beyond the present possibilities. For comparison we have also shown the gluon densities for the GRV 
LO parametrization j^, for the MRS parametrization and also for the pure BFKL evolution as 
calculated in [T^ . 

Putting the found gluon densities into Eqs. (9)-(15) we obtain the (singlet part of) proton structure 
function F 2 {x, Q^). The results are illustrated in Figs. 7 and 8 for the Q^-dependence and Figs. 9 and 
10 for the x dependence. As for the gluon densities, the experimentally investigated region is shown 
separately in Figs. 7 and 9, in which the existing experimental data from ^ are also presented. 



Finally, to see a possible influence of subasymptotic effects, we have repeated the procedure taking 
as a starting point for the evolution a lower value x = 0.001. The resulting gluon distributions and 
structure functions are also presented in the above figures. 

6 Discussion and conclusions 

To discuss the obtained results we have to remember that they involve two quantities of a different 
theoretical status. One is the pomeron wave function (j) which can be identified with the gluon 
distribution (up to a factor) on a rather solid theoretical basis. The other is the quark density (which 
is equivalent to the structure function), for which we actually have no rule for the introduction of 
a running constant and which in the present calculation involves a semi-phenomenological ansatz 
(10). Evidently the results for the latter are much less informative as to the effect of the running 
coupling introduced in our way. Therefore we have to separately discuss our prediction for the gluon 
distribution, on the one hand, and for the structure function, on the other. 

Let us begin with the gluon distribution. Comparing our results with those of GRV, which cor¬ 
respond to the standard DGLAP evolution, we observe that at high enough and low enough x 
our distributions rise with Q and 1/x faster than those of GRV. This difference is, of course, to be 
expected. The hard pomeron theory in any version predicts a power rise of the distribution with 1 jx 
to be compared with (19) for the DGLAP evolution. As to the ^-dependence, the fixed coupling 
(BFKL) hard pomeron model predicts a linear rise, again much stronger than (19). Our running 
coupling model supposedly leads to a somewhat weaker rise. From our results it follows that it is still 
much stronger than for the DGLAP evolution. However one can observe that these features of our 
evolution become clearly visible only at quite high Q and 1/x. For moderate Q < lOGeV/c and/or 
x > 10“"^ the difference between our distributions and those of GRV is insignificant. As to the DGLAP 
evolved MRS parametrization, it gives the gluon distribution which lies systematically below the GRV 
one and, correspondingly, below our values, the difference growing with Q and 1/x. 

We can also compare our gluon distributions with the pure BFKL evolution (fixed coupling) results, 
as presented in [f2]. One should note that the initial values for the evolution chosen in [12] are rather 
different from ours (borrowed from GRV). The initial gluon distribution in [12] is smaller than ours 
by a Q-dependent factor, equal to ~ 2.5 at Q = 2 GeV/c and ~ 1.4 at Q = 30 GeV/c. If one roughly 
takes that into account then from Fig. 3 one concludes that at moderate 1/x our evolution and the 
pure BFKL one lead to quite similar results. However at smaller x (Fig. 4) one observes that our 
running coupling evolution predicts a weaker rise with Qj as expected. 

Passing to the structure functions we observe in Fig. 7 and 9 that our results give a somewhat 
too rapid growth with 1/x in the region 10“^ < x < 10“^ as compared to the experimental data (and 
also to the parametrizations GRV fitted to these data). With the scaling factor 1.3 introduced to fit 
the data at x = 0.01 we overshoot the data at x < 10“^ by ~25%. Without this factor we get a 
very good agreement for x < 10“^ but are below experiment at x = 0.01 by the same order. This 
discrepancy may be attributed either to subasymptotic effects or to a poor quality of our ansatz (10). 
Gomparison with the result obtained with a lower starting point for the evolution x = 0.001 shows 
that subasymptotic effects together with a correct form of coupling to quarks may be the final answer. 
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8 Figure Captions 

Fig. 1 The forward amplitude for a pomeron coupled to a virtual photon. 

Fig. 2 The singlet part of the structure function of the proton at a; = 0.01. The continuous line is the 
result of our calculation while the dashed line correspond to the GRV prediction. 

Fig. 3 The gluon distributions as a function of evolved from x = 0.01 and x = 0.001 for the 
experimentally accessible kinematical region. Standard DGLAP evolved parametrizations (GRV-LO 
and MRS) and the BFKL evolved distributions from [12] (we report only few points connected by 
lines) are shown for comparison. 

Fig. 4 Same as Fig. 3 for asymptotically high values of and 1/a;. 

Fig. 5 The gluon distributions as a function of x evolved from x = 0.01 and x = 0.001 for the 
experimentally accessible kinematical region. Standard DGLAP evolved parametrizations (GRV-LO 
and MRS) and the BFKL evolved distributions from [12] (we report only few points connected by 
lines) are shown for comparison. 

Fig. 6 Same as Fig. 5 for asymptotically high values of and 1/x. 

Fig. 7 dependence of the singlet part of the proton structure function obtained by evolution from 
X = 0.01 and x = 0.001, compared to the GRV prediction and the ZEUS 94 data. 

Fig. 8 Same as Fig. 7 for asymptotically high values of and 1/x. 

Fig. 9 X dependence of the singlet part of the proton structure function obtained by evolution from 
X = 0.01 and x = 0.001, compared to the GRV prediction and the ZEUS 94 data. 

Fig. 10 Same as Fig. 9 for asymptotically high values of and \/x. 
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